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1. Introduction 

The theory of fundamental strings has attracted a lot of interest in the past decade as a 
possible candidate for a unified theory of all the known interactions. Despite heroic efforts, 
crucial questions still remain unanswered. For example, the basic degrees of freedom of 
this theory are still not well understood. In particular near the Planck scale, where we 
expect the stringy behavior to become more evident, it is very poorly understood what is 
the correct configuration space and which states dominate the dynamics of string theory 
at that scale. 

One way to probe this region would be to study a very hot ensemble of strings. There 
is, however, a difficulty as one tries to to heat the ensemble past a limiting temperature 
Th = 13 H^ |T|50 5 the famous Hagedorn temperature. The number of states of string theory 
at a specific energy level increases exponentially with the energy. As a result, the parti- 
tion function is infinite for /3 < fin due to the competition between the entropy and the 
Boltzmann factor e~^^ . What happens when we try to increase the temperature past this 
point is, despite much effort, merely a speculation. 

Since the effective string tension vanishes as one approaches Th one possibility is that 
pumping more energy into the ensemble to increase its temperature merely results in the 
formation of longer and "wigglier" strings by exciting higher oscillation modes [^ . This 
might be related to a "string uncertainty principle" 

he GnE 

^x^ — + -^-^ , 

where as we try to probe smaller distances by scattering more energetic strings, we obtain, 
after a point, poorer resolution due to the formation of longer and longer strings. 

Another suggestion is that the Hagedorn transition is really a first or second order 
phase transition, where above Th we obtain a large genus zero contribution to the free 
energy. In a string theory with smooth world sheets we expect the genus zero contribution 
to the partition function to be zero. The reason is that the partition function is calculated 
by toroidally compactifying the (euclidean) time direction with radius R = P/2n. Since 
the sphere is simply connected, it cannot wrap around the time direction and it should give 
a trivial (3 dependence. This suggests that smooth Riemann surfaces are not appropriate 
for describing the high temperature phase of the string. The following picture is suggested 
to hold [^|-§1: The modes f and f* that wind once around S^ become tachyonic at 
Th and are supposed to acquire nonzero expectation values < (p > and < (f* > around 
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which we should look for stable solutions (in fact due to the dilaton the transition occurs 
at Tcrit < Th). Then the (p vertex operator creates tiny holes on the world sheet that 
wind around the time direction tearing the world sheet. This suggests that the degrees of 
freedom in the high-T phase are drastically reduced. In fact the large T behavior of the free 
energy implies that there are fewer fundamental gauge invariant degrees of freedom than 
any known field theory that lives in the embedding space (for example for the closed strings 
one finds that the number of degrees of freedom is the same as that of a two dimensional 
field theory!). 

At this point it is useful to consider the parallel with low energy QCD. It has been 
claimed long ago that large A^ QCD or the strong coupling limit of QCD p can be 
formulated as a string theory (for a recent review see and references therein). In favor 
of this point of view is the early success of the dual resonance models in describing the low 
energy meson resonances. There are, however, difficiilties with this picture. For example 
the string in the relevant non-critical dimensions is not Lorentz invariant in the light 
cone quantization and has one too many oscillators because of the Liouville mode in the 
Polyakov quantization. 

It is believed that SU{N) QCD undergoes a first order deconfining transition at 
some temperature Tdec- This can be the analogue of the Hagedorn transition of string 
theory. We understand the high energy regime because of asymptotic freedom. The basic 
degrees of freedom are weekly interacting quarks and gluons. In the strong coupling limit 
perturbative QCD is no longer a good description anymore, since the perturbation series, 
even as an asymptotic series, is not a good approximation. Hadrons, mesons and glueballs 
are the basic degrees of freedom and an effective string theory description in terms of 
thin chromo-electric flux tubes might be possible. The low temperature phase is described 
by smooth world sheets and it is expected that such a description breaks down in the 
high temperature regime; it is the underlying field theory that gives the basic degrees of 
freedom in that case. In string theory, however, a 26 or 10 dimensional field theory would 
be inappropriate for the high temperature phase, since such a theory would not have a 
good ultraviolet behavior. 

The study of finite temperature strings also suggests some perhaps deep connection 
between multiplicative number theory and statistical systems with an exponentially grow- 
ing density of states [p!0| -p!3|. There exist several examples, like the bosonic, fermionic 



or parafermionic Riemann gas and the known string and conformal field theories, where 
the partition function is related to well studied number theoretic multiplicative arithmetic 



functions and modular forms. Understanding such a connection for physically interesting 
string theories could be of great importance for calculating and studying the behavior of 
relevant physical quantities. The number theoretic methods developed are quite power- 
ful and can be used to calculate the partition function or give a very good estimate of 
the asymptotic behavior of the level densities of the above models. We should also not 
underestimate the possible mathematical importance of these questions. 

This presentation is organized as follows. In section two we introduce the main rele- 
vant concepts. In section three we present the connection between string theory at finite 
temperature and number theory. We stress the relation between the level density p{N) 
of the bosonic string, the Dedekind eta function ?7(r) and Ramanujan's r-function r(?i). 
We explain how one can use Hardy and Ramanujan's results to compute p{N). In section 
four we make this connection deeper by discussing the Riemann gas of bosons, fermions 
or parafermions. In section five we discuss the micro canonical and canonical ensemble 
of strings emphasizing the limit of validity of the canonical ensemble near the Hagedorn 
temperature. In section six we present the derivation of the asymptotic density of states 
for the micro canonical ensemble and the effect of introducing chemical potentials in order 
to impose conservation laws in the computation of Ph- In section seven we describe how 
we can obtain the partition function with the string path integral on the torus. In sec- 
tion eight we discuss the possible physical behavior of the string ensemble near Ph and in 
section nine we discuss the results of Kogan and Sathiapalan and Atick and Witten. 

2. Strings as a Statistical Mechanical System 

Consider a d dimensional string theory described by the embedding 

of the Riemann surface (complex curve) E with coordinates (a, r) into the d dimensional 
Minkowski space R^~^'^ with cartesian coordinates X^, fx = l,2,...,(i. The classical 
free equations of motion 

{dl-d'^)X^ = (1) 

and the boundary conditions, which for the open string are X'^{a,T) = for a = and 
TT, determine the mode expansion 

x^{a,T) =a:(^+p^r + z^-^e-'"^cosna. (2) 



Upon quantization, the constraint (Lq — a)\x >= on the physical states determines the 
spectrum of the theory 

-M'^ = N-a. (3) 

In the above formulas Lq is the zero mode of the stress energy tensor, a a normal ordering 
constant and N = Yl'^=i na-n ■ otn = Xl^i ^<^n ' <^n is the level number operator. The 
number of states at a particular level increases very rapidly with A^ and as we will soon 
see, it increases exponentially with A^ for large N. In order to compute the number of 
states at level A^, we consider the generating function for the level degeneracies 

F{z) = Tr ^^ 

N=0 

where d{N) is the number of mass eigenstates at level TV. We can compute F{z) using 



X]*^""' 



F{z)=Ttz^ = IJTr 2''""-"- 






|d-2 



d-2 (5) 



where ^ f{z) = Yl rrpr = Yl p{n)z'^ is the classical partition function of Euler. It is the 
generating function for the number of unrestricted partitions of A^ into positive integers. 
Proof: Since 

oo 

f{z) = l[{l + z^ + z^^ + z'^ + ...), (6) 

r=l 

a typical term 2;^ in f(z) is obtained by taking one contribution from each factor r = 
1,..., 00 such that 2^ = (2)^i(z2)^2(^3^fc3___(^n^fc„ withk = ki + 2k2 + 3ks + ...+nkn. This 
yields a partition of k into sets {ki}. 

3. Singularities of f{z) and Asymptotics of p{N) 

The function f{z) has an essential singularity at all rational points of the unit circle 

2TTip 

S i.e. the set of points {e 1 : p^q E Z} = {zp^g}. The essential singularity arises 



' The tachyon may be accounted for by modifying this to F{z) = -|/(2;)| '^. 
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because Zp=l for an infinite number of integers n = /c ■ gV/c G Z. In order to estimate 
the asymptotic behavior of the density of states we have to study the behavior of f{z) as 
z -^ 1~ . A crude estimate gives 



oo ^ / oo 

hm f{z) = hm I I = hm exp — ^ ^ ln(l 



zn 



Ti=l \ n=l 



- £m exp [J: ^^ = £m exp J^ . (7) 



v^?Ti,n=l 



— ', m(l — 2^^ 

1=1 ^ 



oo 



hm exp y — 

z^i- \ ^^-^ m^ 1 — 

1,771=1 



This is the first appearance of exponential growth in the problem. One can improve the 
above estimate by using the modular properties of f{z). Writing z as e'^^*'^, the classical 
partition function takes the form 

oo ^ 

/(^) = n 13^^ ' (8) 

n=l 

which is the partition function of a single boson on the torus with modular parameter r. 
This is the first connection we see with additive number theory (combinatorics) and 
conformal field theory. The function f{z) is closely related to an important modular form, 
the Dedekind eta function ?7(r) via 



oo 



rj{r) = e^ n (1 - ^'""1 = ^*//(^) • (9) 



n=l 

This function has many fascinating number theoretic properties. For example 

24 _ 



oo 



v'\r) = j^^=z{l[{l-z-)\ =$:r(n).- (10) 



, Ti=l J n=l 



defines Ramanujan's r-function T{n). 

We list below some interesting properties of r(n) 
1) Ramanujan conjectured that T{n) is a multiplicative function, i.e. 



T 



(mn) = T{m)T{n) if (m, n) = 1 



where (m, n) is the greatest common divisor of m and n. This was proved by Mordell 



2) T{m)T{n) = E d^^r{^). 

d\ (m,n) 

3) If t{p) = (mod p) then r(pn) = (mod p) Vn. 

4) Dyson []TB| proved that 

, , ^^ {a — b){a — c){a — d){a — e){b — c){b — d){b — e){c — d){c — e){d — e) 

'^^''^ ^ ^ 1!2!3!4! 

a,b,c,d,e 

where a, 6, c, d, e satisfy 

a, 6, c, (i, e = 1, 2, 3, 4, 5 (mod 5) , 

a+b+c+d+e =0, 

and 

a^ + 6^ + c^ + rf2 + e^ = lOn . 

The function rj{r) is a modular form of weight +| since 

r7(-l/r) = (-zr)i/2r?(r). (11) 

Using ( [TT] ) we can map the region 2; — > 1~, where /(z) has singularities, to the region 
z' — > 0"*", where f{z) is regular and close to 1. Under r -^ — I/t, z = e^'^*'^ ^ 2' = e"^'^*''^ 
and since z' = exp ( j^ ) , we obtain 

Thus 

lim /(2;) = lim ( ^-^ ) z^ exp ( — ^ ) (12) 

Note that (p!2D is a slightly improved estimate of the behavior of f{z) as z — »• 1~ as 



compared to (^. 



To extract p{N) from the above results we follow Hardy and Ramanujan |jT^ who first 
applied the techniques of complex analysis to the combinatoric problem of determining the 
asymptotic behavior of the number of unrestricted partitions of an integer A^. 



where the contour is taken around the origin and within and close to the unit circle. Then 



1 /" , / — Inz 



1/2 



^' _(Ar + ^)ln^ 



61n2 ' 24' 

A saddle point evaluation of p{N) (see appendix) leads to the result 



(14) 



A better result for p{N) was derived by Hardy and Ramanujan. They showed that 

^^^^' 2nV2dN\ Xj 

(16) 
V3 /2,^ 1 \ d /e3=^^\ 

H 7= cos -Nn n -— —^ + • ■ • 

nV2 V3 18 J dN ^ ' 




where Xn = \/N — ^ and c = ttW |. For A^ = 100 and 200, summing 6 terms of this 
series gives p(lOO) = 190,569,291.996 and p(200) = 3,972,999,029,338.004, whereas the 
exact results are p(lOO) = 190, 569, 292 and p(200) = 3, 972, 999, 029, 388 respectively [|l|]. 



Later an exact analytical result was found by Rademacher [|T^ using a slight variant of the 
Hardy-Ramanujan analysis 

r'W-^,Y.MN),i^(-^y (17) 

where Aq{N) is a certain sum of roots of unity. 

4. Quantum Statistical Mechanics and Multiplicative Number Theory 

Multiplicative number theory in the context of statistical systems with exponentially 
growing number of states has been considered first by Julia and Spector[lTO|-[TT[|. Consider, 
instead of the usual Fock space of a quantum field theory, a different labelling of the 
quantum states. Place the creation operators {a]} in cardinal order and associate the 
^th pj-jjjie number with the i^^ creation operator. Then label the state rii('^D'^'|0 > t)y 
the integer A^ = 11^^?- Since the factorization of A^ into prime numbers is unique, this 
provides a basis for the Hilbert space of states. This is known as the Godel numbering 
|T^ . We are arithmetizing this way a free quantum field theory the same way as the formal 



calculus of propositions is arithmetized by associating the i^^ symbol in a proposition with 
the i*^ prime number and determining the corresponding multiplicity by the symbol itself. 
A simple example of a system with exponential growth of the number of states is the 
Riemann gas. The partition function of such a system is defined as 

Z{l3)=Tve-^^, (18) 

with energies Ei = In pi. Then 

Zi/3) = ^ exp(-/3 ^ Ti In pi) 

oc (19) 

N=l 

is the (ordinary) Riemann zeta function. The partition function therefore, has a simple 
pole at /3 = 1. This can alternatively be seen by computing the density of states 

p{E) = -l^p{N)A{N) = 

-i^ky-^(i + ^ + ---)- (20) 

= exp(£;) + ^ + 0{e-^) . 

We see that (3h = 1 corresponds to the simple pole of ({(3) at /3 = 1. 

Generally, multiplicative generating functions of the form (Dirichlet series) 

FW = ± i^ (21) 

n=l 

are of interest, where f{n) is a general multiplicative arithmetic function on the natural 
numbers valued in some field (usually the field of real or complex numbers) . As in the case 
of the ordinary Riemann zeta function, F{P) generally has a simple set of singularities but 
a rich structure of complex zeroes. 

The case f{n) = ii{n), where ^(n) is the Mobius function, gives the partition function 
of the fermionic analog of (^8]). This gives Zp{(3) = ^^^- Similarly the system of k- 
parafermions is associated with Jir 



/(n) = ^fc(n) = ^ 1 if^ = n.bO'^^withO<r,</c-l; ^22) 

otherwise. 



It gives 

Parafermions of order 2 and cxd correspond to fermions and bosons respectively. Note that 
a representation of the ordinary Riemann zeta function can be obtained by tensoring an 
infinite set of parafermionic gases of order k at successively lower temperatures, such that 

CO 

aP)^Z^{P) = l[Zk{k'^P). (24) 

n=0 

We hope that the above exposition has given a flavor of the important connection between 
arithmetic gases and multiplicative number theory. Much remains to be understood for 
the case of statistical systems that arise from string theories and conformal fleld theories. 

5. The String Microcanonical and Canonical Ensemble 

We are now ready to compute the microcanonical density of states for the string 
ensemble. Since the mass spectrum is given by 

-M^ = N - const , (25) 

the asymptotic form for p{N) translates into the string density of states 

0(M) ~ M"" exp (/3M) . (26) 

The microcanonical density of states Q{E) is deflned by 

n{E) = Tr 6{E- H) 

^^Tr e'^^^~^^ 



2n 

/3o+ioo 



^e^^Tr e-^^ (27) 

2/TZ 



Pq—ico 
/3o+ico 

d(3 



2m 

where (3 = i(3 and (3o is chosen for the convergence of Z{(3) = Tr e~^^, the canonical 
partition function. Inverting the Laplace transform we obtain 

/•CXD 

Z{(3) = / dEn{E)e-^^ (28) 

Jo 
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In a typical statistical mechanical system with A^ degrees of freedom Q{E) ~ E^ for large 
E. Thus the entropy S = lnO(i?) ~ NlnE is extensive. For large A^, the integrand is 
sharply peaked and the canonical partition function can be evaluated by a saddle point 
approximation. Since 

/•OO 

Z{(3)= dE exp {In n{E)-pE) (29) 

Jo 

the saddle point occurs at Eq such that 



dlnO(E) 



dE 



E=E, 







Eq 



(30) 



i.e. Eq = -Q = N T. Thus fixing the temperature T (canonical ensemble) is equivalent 
in the thermodynamic limit to fixing E in the micro canonical ensemble at NT and T is 
really the average energy per degree of freedom. The fluctuations are given by 

OO 

Z{P) = ZoJdE exp 1^{E - Eof-^iS{E) - PE) 



Zo fdE exp!^^iE-Eo) 





' d'S{E) \ 

dE^ J ■ 



(31) 



This may be rewritten in terms of the micro canonical (formal) specific heat Cy- Since 

Sm = lnQ{E) and Tm = (f^) ? we have 



dE^ 



1 dT 
T^dE 



1 



1 



V 



T2 {Cv)m 



(32) 



where Cy = ^. Then 



[Cv] 



1 /d^S 



M 



T2 \dE^ 



(33) 



Thus (^rp becomes 



Z{(3) = Zoexp{-^{E~Eo)':^{C- 



v)m 



(34) 



Convergence requires that Cm > in order that saddle point approximation is valid. 

For strings Q{E) ~ E~'^ exp PmE and as one approaches Th from below the integrand 
Q{E)E-f^^^ fiattens out. Typically a > 1 and Z{(3) is weU defined for (3 = (3h, but it 
diverges for (3 > 13 h- Fixing T as /3 ^ (3^ no longer corresponds to fixing a precise E and 
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the canonical and micro canonical ensemble need not agree. The entropy S = In fl{E) 
(3hE — alnE is no longer extensive as well. For further discussion see section 7. 
For free particles 

Z{I3) = n (1 - e-^^^i'')"' n (1 + e-^^^1/) . 

kib kif 



Therefore 



oo ^ 

lnZi(3)=J2l 

r=l ^ 



MM _ ^_^yMM ) (35) 



where the single-string partition functions /b and fp are 

oo 

k,b Q 

OO 

fp = 5^e-^^-'/ ^ fdE e-^^LOFiE) 



(36) 



Substituting (|3^) into (|23) we obtain the multi-string density of states from the single 
string density of states 



^" ' oo 



.,.,^ J ^e-exph:(MM-(-:rMM) . (37) 

The r = 1 term gives the familiar Maxwell-Boltzmann expression 

oo ^ « n 

n{E) = 5^ - / n ^^^ u;{E,)S{E - J2 Ei) (38) 

where a;(£') = ujb{E) + ujf{E) 

6. Prom the One-String Density of States to the Multi-String 
Density of States for Toroidal Compactification 

Consider d-dimensional strings compactified on an internal manifold of dimension c 
that is for simplicity taken to be a torus with radii Ri with i = 1, 2, . . . , c. In computing the 



full density of states, it is important to incorporate the effects of conservation laws [19| such 



as conservation of momentum, winding number and, in the case of the heterotic string, of 
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U{1) charges or windings/momenta in the internal directions. For every conserved charge 
Qa we introduce a chemical potential fiA and compute 

Q{E, ixa) = tr le^'^'^^Q^SiE - H)] 

(39) 

The condition that the total conserved charge be some fixed value qa (typically qa = 0) 
is then enforced by multiplying Q{E, ^a) or Z{(3, ha) by exp{— 27rZ|Uyiq'yi} and integrating 
Ha over the range (— |, \) 

Vt{E,qA) = j\ diiA^{E,HA)e-^^'''^''^ ■ (40) 

The full density of states will be derived from the single-string density of states as discussed 
in the previous section. In order to compute the single string density of states, consider 
the energy of a general string state 



*^+i:( 

i=i ^ 



2 



E =k' + } A^ + i'^R^r^^] + 4nL + 4nfl . (41) 

Here k is the momentum in the non-compact directions and rii and rrii are the momentum 
and winding quantum numbers, respectively. The integers ul and ur label the internal 
level numbers describing oscillatory mode excitations of the string. 

As described in section 1, the number of left-moving (right-moving) excitations at a 
given level is given by the generating function F{z) = Y2^ ^'^La (^F{z) = ^^ ^"^««) where 
the sum runs over all possible states a. For a bosonic sector f{z) = rj~'^^{z) = -/{z)"^^ (the 
extra factor ^ corresponds to the tachyon), while for a superstring sector f{z) = ^^\ . 
As before we compute the asymptotic level densities using a saddle point evaluation 

d{n) = -i. / -%f{z) . (42) 

The resulting level densities for left and right sectors are 

d{nL) ~ n;^^^+'^/''e4"v^^^^^ , d{nR) ~ n^^^'^+^^^^e^''^^^^^ (43) 

where aL,R = ^ ior a. superstring sector and aL^/j = 1 for a bosonic string sector (so for the 
heterotic string a^ = 1 and gr = |). Physical states must also obey the level- matching 
condition Lq — Lo = 0, or ul — ur = — J2i ^i^i- This and (^l]) give 

2 



Snr = E^ - k"^ - { ^ + 2Rim 



n 
iL=E--k--{ 

. 2 

SuR = E'^ -k'^ - { — - 2Rim, 
Ri 
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Therefore the single-string density of states is (including the chemical potentials k, h and 



z^ 






mi.Tii 



(44) 

I 4:TTy/ a LnLiE,k,ni,mi)+4:iTy/aiinii(E,k,ni, mi) \ 

riLiE, k, 7i„ m,)(^+i)/4ni?(£;, /c, n,, m,)(<^+i)/4 • 
Here L> = (i — c— 1 is the number of non-compact space directions, and Vd is the volume 
of the non-compact directions.) 

Note that for each Ri there are regions with different behaviour. liRi^ E ^ ^/Ri we 
ignore the winding quantum numbers rrii and the sum over the closely spaced momentum 
levels is well approximated by an integral over a continuous spectrum of momentum; thus 
the string behaves as if the direction i were noncompact. Likewise, if 1/Ri ^ E ^ Ri the 
sum over momenta drops out and the sum over winding becomes continuous; this situation 
is dual to the previous one under R -^ 1/2R. Finally we have the truly high-energy regime, 
E ^ Ri, l/Ri- Here we can closely approximate the sum over m^, rii by an integral drriidni. 
Furthermore, for large E we can evaluate the integral (|3^) by saddle point methods. If 
we expand the arguments of the square roots to first order in /c^, {rii/Ri ± 2Rmi)^ (valid 
for large E in the saddle point approximation since the region {ui/Ri ± 2Rmi)'^ / E'^ ^ 1, 
k'^/E'^ ^ 1 dominates), the integral becomes quadratic and yields 

u;{E, K, ^, v) ~ -^e^-^'^-'^'-)^ (45) 

where 



.2 , „ „Rx2 



( ^2 (j^s]^Y (j^^eRV 

13h{i^, /U, J^) = Ph - V2n < — — -= -I -= \ -= — 

/aL + \/aR Jal Wor 



(46) 



and 

Ph = V2n {^/aZ -f a/or) 



is the usual Hagedorn temperature 0. Thus in effect the Hagedorn temperature depends 
on the chemical potentials. 

This result can now be used to derive the multi-string density of states as sketched in 
the previous section. First we will compute the Maxwell-Boltzmann contribution (|3^) to 
the total density of states, and then we will argue that the corrections due to Bose-Einstein 
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or Fermi-Dirac statistics do not substantially alter the asymptotic form of the density of 
states. 

The MB expression ( |38| ) becomes a sum of terms of the form 

^/nTw^''"'"''"""'«(^-Es.i ■ (47) 

i=l i \ i 

The integrals over Ei are divergent at the lower end and must be cut off at some energy 
mo>Ms where the asymptotic expression ( ^Sf ) is no longer valid. It is clear that the 



dominant contribution to (|47| ) is that in which most of the energy is in one string, Ei ^ E. 



Thus the integral is approximately 



1 p|3H{|^,^J.,y)E 

- [ln(£;/mo)]""^ for L> = 0. 



(48) 



(n-1)! E 

Summing over n then gives 



n{E,^,^,jy)r^ ^^^^^^ e^lV''^- iforl^>0; 
mo 



(49) 



The density of states without the conservation constraints is given by setting k = \x = 
z^ = 0; this expression has been found in the literature. To obtain the density of states 
including the constraint of zero total momentum and winding, we use (|i6|) and integrate 

over K, (U, V to find 

pPhE 

where we define ^d = 1 if D = and ^i) = otherwise. Imposing the constraints changes 
the power of E\ this is in contrast to other thermodynamic situations where conservation 
laws are not so important essentially because one string carries most of the energy and 
bears the burden of the conservation laws. 

We can also compute corrections to ( |50D due to Bose-Einstein and Fermi-Dirac statis- 
tics. One finds that the effect of these corrections is to multiply (|50|) by a function that 
depends on low-energy physics but to leading order not on E. 
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Thus for the type II string the asymptotic density of states is 

^(E) - ^.^ . (51) 

If we consider the heterotic E^ x Eg or Spin(32)/Z2 string we see that conservation of 
the 16 U(l) charges should also be imposed. In the bosonic formulation this is equivalent 
to momentum/winding conservation for the internal degrees of freedom, and, as one can 
easily convince oneself by a slight modification of the above argument, gives an extra factor 
of 1/E^. Thus the asymptotic density of states for the heterotic string is 

^(E) - ^,33^ . (52) 

It will be useful for us to deduce which string configurations contribute most signifi- 
cantly to the total string density of states; this can be done by inspection of (|47| ) or the 
corresponding expression including quantum statistics. As we have already noted, the 
dominant configuration for the integral for a fixed number n of strings is when Ei ^ E, 
i.e. when most of the energy is in a single string. For D > it also appears that string 
configurations with a small total number of strings dominate the sum, but for D = 
configurations with a large total number of strings dominate the sum. (These qualitative 
features also hold true when corrections due to quantum statistics are included.) We will 
return to a discussion of this single string dominance of the energy in section 7 

7. Path-Integral Derivation 

In the meantime, we turn to the second derivation of the string density of states 
pOH2^ and 0. Once again {^7^ is used to relate Q{E) to the complex function Z{(3). 
Next recall that the free energy F = — 4 In Z for a single free particle of mass M can be 
computed from a first-quantized one-loop path integral on a space with compactified time 
direction of circumference (3. The action is Sp = ^ J^ dT[e~^ {dX / dr^ + eM^], and in 
the path integral we integrate over all maps X^{t) of the circle into the target spacetime 
S^ X R'^~^ and over all one-metrics g = e^ on the circle. Explicitly, 



n= — CX3 



where in the path integral over metrics we must eliminate the volume of the diffeomorphism 
group, and where the boundary condition on the X path integral is -'^^(l) = X^{Qi) + 
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n/3, X'^{1) = X*(0); n is the number of times that the circle winds around the compact 
time direction. We also define F to be zero for bosons and one for fermions; the extra 
factor (—1)"^'^ corresponds to anti-periodic boundary conditions for fermions in the second- 
quantized formalism. One easily shows (taking into account the Killing vector on the circle) 

that 

Vg p ds 



Vol(Diff) Jo s 
where s is the proper length of the circle. One can also show 



J n 



2 



where V is the spatial volume. These give the free energy 



^ri?£(2..)-/^ f e-4|i-*M)-. (54) 



2 JO 

^ n= — oo 



It is fairly simple to demonstrate that this is equivalent to the standard expressions for 
the free energy of a single free particle. 

Likewise, the free energy for string is obtained by doing the functional integral over 
maps of the torus into 5"^ x R'^~^, with the (conformal gauge- fixed) action which we may 
take to be that of the heterotic string (our units are such that a' = ^) 

S=l-Jd^a {d^X^'d.X^ + rd,^^ + r^-.^^. + A^a.A*) . 

(The \i provide a fermionic representation of the gauge degrees of freedom.) The maps 
are allowed to wind in the time direction on the torus, but not in the space direction, as 
in the case of the particle. The result is 



V _5 r dT2 p ^ ,^ ._5 1 

— 71 ^ / dTi(2nT2) ^ , 

4 Jo r2 y_i viQ^viq)' 



F^--^ I —I dT,{27vr2)-''———-j^ch[G]{q) 



■Kl,TT2=0,l n ^ ^ ■' 



(55) 



/j^NnTTl 



Here q = exp{27rz(ri + zt2)}, t] is the Dedekind function; c/i[G] is the character function of 
the gauge group G = Eg xEg or SO(32), ch[E8 x Eg] = {^l+'&l+'&l)^ /Atj^^, ch[SO(32)] = 
{■dl^+'&l^+'df)/2r]^^; and Cqo = ^|, 6*01 = -^j, Cio = -??|, and Cu = 0. This expression 
is analogous to (^4]); the integral over nT2 corresponds to that over s, the sum over n is the 
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sum over sectors where the time direction on the torus wraps n times around the target 
time dimension, and the rest of the expression integrated over ri gives the sum over aU 
single string states of (— l)"^exp{— M^7rr2/2}. We could therefore alternatively derive 
this expression directly from (Q) by thinking of the string as a collection of particles 
corresponding to its various modes. 

To make this expression look more familiar from the string point of view, we recall 
that the region S of the r = ri + iT2 plane |ri| < |, < r2 < oo corresponds to an 
infinite number of copies of the fundamental domain JF defined by |ri| < |, |r| > 1; S 
is tiled by images of JF under maps r' = {pr + (l)/{rr + s), p, q, r, s G Z, ps — qr = 1. 
The free energy for the heterotic string can therefore be rewritten as an integral over the 
fundamental domain 






\m7r2+n7ri +77171 



^expl -^^ ("^^T-| + {n-mrif) \ (-1)' 

m,n ^ ' 

and likewise for the type II string 

nnrt nn ^ ^ 



(56) 



— TV ^ 



(57) 



m,n 



From ( |56| ) one can proceed to calculate the asymptotic density of states . 

In these expressions the sum over m is the sum over maps of the torus in which the 
space direction of the torus is wound m times about the target 5"^ . While one might think 
of these as corresponding to strange states of the string in which it wraps around the time 
direction, this interpretation is not necessarily correct if the fundamental expression for 
the free energy is (^5[) which directly corresponds to Z((3) = tre~^^. This trace contains 
only ordinary string states and not the strange 'winding states;' the latter appear to arise 
only as mathematical artifacts when we transform the expression for InZ as an integral 
over S into an integral over JF. 

Note the following interesting point. (^) is equivalent to the modular invariant ( [5^ ) 
the cosmological constant A = 0. However, if A 7^ 0, (|55|) gives 00 ■ A + F ' and (|5^) gives 
A + F ' where A is cosmological constant as it is usually defined in string theory (as an 
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integral over JF) and F ' is the temperature dependent part of F. This is also a difficulty 
for string field theory -i.e. to naturally produce the fundamental domain. 

8. Physical Behaviour of Strings at High Temperature 

We now turn to an attempt to extract a better understanding of the physics of strings 
at high temperature from the above calculations 0,^- To start with, assume that we 



have at our disposal a thermal bath that can interact with the string ensemble to heat it 
up. Within the context of fundamental strings this assumption is not necessarily valid as 
all matter is presumably made out of strings, but we make this initial assumption to aid 
us in isolating the source of any strange behavior that we may encounter; we will see that 
the assumption of the existence of such a reservoir doesn't obviously appear to be related 
to any such behavior. Therefore suppose that we have a heat reservoir with a density of 
states Qr{E). 

The total density of states for the string ensemble in contact with this reservoir is 

VLt{E)= f dEsVL{Es)VLr{E-Es) . 
Jo 

If we assume that the reservoir is large (in a sense to be made precise shortly) , then most 
of the energy will be in the reservoir, i.e. E ^ Eg. In this case we can expand 

lnO,(E - Es) = InQriE) - Es-^\nQr{E) + ^£;2^ lnO,(E) + - • • • . 

We identify -S^ \\iVtr{E) = /3 as the inverse temperature of the reservoir. We can drop the 
higher terms when the reservoir is sufficiently large; precisely, E'^-^jSriE) <^ 1 or 

l.Cy.^y^»E^.. (58) 



Thus we can use the ensemble 



/"OO 

rtriE) = QriE) / dEse~^^'n{Es) 
Jo 



(59) 



as long as we ensure that the heat capacity of the reservoir is large enough in the sense of 
; replacing the upper limit by oo has no significant effect. 
Now we are prepared to study the properties of the string gas. As we raise the 



temperature, (^9|) stays well defined all the way to (3h- Even at (3h (H^) is well defined 
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because of the power law tail ~ 1/E^'^~^^' . (Of course Vt{E) takes a different, convergent, 
form for low energy; since the string is a gas of massless particles at very low energies, we 
expect VL{E) ~ exp < a^^E^'^^^^'VjJ \ there for some constant a). Furthermore, 

the mean energy density 

is finite for (3 = (3h', e~^"^Vt{E) will presumably look something like Fig. 1 where the 
maximum Eq is somewhere between the low energy and high energy domains where Vt{E) 
is known, i.e. £q will be of order one in string scale units. 

What happens when we attempt to raise the temperature further? It is apparent 
from (|59D that for (3 < (3h this ensemble is not defined; its mean energy becomes infinite. 
Therefore it seems that we cannot pass Ph- But in principle we could imagine adding more 
energy to the string system in an attempt to raise its temperature; where does this energy 
go? A clue was pointed out at the end of section five: at large energies the density of 
states is dominated by one (or a small number of) highly energetic, and therefore long and 
wiggly, strings. Thus it seems quite reasonable that any energy that we add to the system 
in an attempt to raise the temperature goes into making strings in very high oscillation 
modes. The low mass modes are populated with a thermal distribution at /?//, and the rest 
of the energy is in oscillation modes. It may be useful to think of this transition to long- 
string dominance as a higher-order phase transition with an infinite latent heat. The low 
temperature phase is what we know, and the 'high temperature phase' is the configuration 
where the fractional energy in long string is one, but this phase is inaccessible because it 
corresponds to infinite energy. (This is assuming that the fundamental degrees of freedom 
are not drastically modified at some finite energy density, as could well be the case.) 

9. X-Y Model and the Berezinskii-Kosterlitz-Thouless Transition 

Consider the string action in the conformal gauge 

S=— [ dadr {dc^X'^d^X'') , (61) 

27r J 

where we take a' = 1/2. Compactifying in the time direction, in order to account for the 
finite temperature, and using periodic boundary conditions for closed string X^((j+7r, r) = 
X^(a", r) and X^{a, t + tv) = X^(a, r) + 7i/35q , we obtain the mode expansion for the time 
coordinate 

X\a, r)=X^ + ^ -f !^ + (oscillators) (62) 

P TV 
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and 

_£^ = iV. + -(- + ^j -l+iV„ + -(^_-^j -1 (63) 

with the constraint Nl — Nfi = m ■ n. The state Nl = Nn = with n = and m = ±1 
has 

£;(T)2 = M(T)2 = -8 + -l;^ (64) 

For T < Th = l/27r-\/2 its mass-squared is positive but at T = Th it becomes massless and 
for T > Th tachyonic. CaU the m = ±1 winding states ip and ip* respectively. Thus the 
Hagedorn divergence (T > Th) corresponds to the existence of a tachyon and a divergence 
at the boundary r2 ^ oo of the moduh space in the path-integral language | 



F{(3) ~ {const.) / dT2 e-^'^^)^2 -^ oo for M'^{(3) < . (65) 

The action for the ip mode given by X^ = ^p> with (^^(a, r + tt) = (/j((j, r) + 27r is 

The action Sip is that of the Villain model which is the continuum limit of that of the X-Y 

model (planar magnet) 

S = - 
2 



5 = -- ^ cos((^i - (^i+5) ■ (67) 



<i,5> 

The vortex solutions of the model are given by 



•p^z) = m In 



z - Zq 
z - Zo 



(68) 



where zq is the position of the vortex and m is the winding number. The corresponding 
free energy is 



-— m - 1) In — 
2 \a- 



Fm = (^m^ - 1) In ( - ) , (69) 



where a is the lattice spacing cutoff. The above model has an infinite order Berezin- 
skii-Kosterlitz-Thouless phase transition aX Qc = — or (3c = 2v27r = /3//! In the low 
temperature regime [g > Qc) the system is dominated by spin wave excitations and is in an 
unscreened scale invariant phase with power law spin-spin correlations and tightly bound 
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vortex-antivortex pairs. In the high temperature regime {g < Qc) the topological order is 
destroyed, the vortices unbind and are free (vortex condensation). The phase is screened 
with exponentially decaying spin-spin correlations characterized by a coherence length ^ 
which breaks the scale invariance. 

Consider a chemical potential fi for vortices, where /x corresponds to the energy of 
dissociation of a vortex-antivortex bound state. The RG flows for the system are in the 
fugacity y {y = e~^'^'^) and T plane as shown in Fig.l. 

For T < Tc and y small the system flows to y = or to /U = oo and vortices are bound 
(region I). For T > Tc and y small the system flows off to y large out of the domain of 
perturbation theory but ^ is effectively small and vortices can unbind. 

To describe string theory in this phase one would have to look at the coupled (3- 
functions for all the massless moduli at /3 = (3h which include g^u^B^^, the dilaton a 
and ip. This might have solutions corresponding to a conformal flxed point for T > Th- 
Atick and Witten suggested that ip gets a v.e.v. < (p > before Th thus obviating 
the difficulties of the Hagedorn divergence. But they found the effective potential for ip 
obtained from its interactions with the dilaton a is of the type 

Veffip) ~ m^{T)p*p - X{p*pf' 

with 

A>0 and m2(T) = -8 + I/tt^T^ . 

This implies a first-order phase transition. Vf,ff{p) must, however, be stabilized by higher- 
order terms and < p > is not calculable in perturbation theory. 

Atick and Witten argued that, had A been negative, there would be a second order 
phase transition at T = Th- For T just greater than Th then 

I 2 I 

,2 \m \ 
\<p>>\ = ^—^ 
2A 

and the free energy is 



F = V{<p>>) 



1 2 

-TO —ra 



4A 4A^2r 

This is a genus zero (sphere) contribution to the free energy and implies the world-sheet is 
spontaneously tearing to become non-simply-connected. They then argue that this picture 
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is valid even for the first-order transition case A > 0. This is not, so far, justified by any 
calculation and fails in the weak coupling limit (7 — > 0. 

Appendix 

In this appendix we compute the asymptotic form of the level density p{N) which led 
to ([T5|). From (^ID we have that 

1 /" / 1 \ ^/^ 

with g{z) = ^ - {N + c)lnz and c = 23/24. Then g'{z) = ^^^ " i - ^. which has 
a maximum at In 2:0 = — 1= ■ —r==- Then 
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g{zo) = -j=VN + c. (71) 



Thus the leading behavior of p{N) is 



p(iV)~exp|.,/|yiV+|^ P2) 



~ exp (Pc^) 
with /3c = ttW |. Corrections are obtained by looking at the quadratic fluctuations 



g{z) = g{zo) + liz- ^o) V(^o) + • • • , (73) 



where g"{zo) = ^^{N + c)"^/^exp (ttw |^== j. Then ( |70|) is a Gaussian integral which 



gives 



''''^' = 37!^''"''rV5^'^l- <"" 



For the string in d dimensions ( |7^) easily generalizes to 



p(N) = (const.) N-^'^+^^/^ exp | nJ^{d-2)VN ] . (75) 
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Figure Captions 

Fig. 1. R.G. flows for the Villain model. The x axis is Tc/T — 1 and the y axis the 
fugacity e~^'^'^ . The x axis is a line of trivial fixed points that can be reached 
from region I. Points from region II fiow to large y and out of the domain of 
validity of perturbation theory. 
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